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ABSTRACT 

The  power  spectral  density  of  a  pulse -position  modulated  (PPM) 
waveform  is  determined  in  part  by  certain  statistics  of  the  modulating 
process.     In  this  paper  a  mathematical  model  is  developed  to  determine 
the  effects  of  three  different  modulating  processes  upon  the  power  spec- 
tral density.      The  processes  considered  are  the  Gaussian,    the  exponen- 
tial,   and  the  uniform  processes,   with  statistical  independence  assumed 
between  the  sample  values  that  modulate  the  PPM  waveform.      The 
Gaussian  process  is  also  developed  for  the  case  in  which  statistical  in- 
dependence may  not  be  assumed;  that  is,    a  correlation  exists  between 
the  sample  values  of  the  modulating  process.      The  model  will  also 
determine  the  power  spectral  density  of  pulse-amplitude  modulated 
(PAM)  waveforms,   or  for  a  combination  of  PAM  and  PPM. 

The  mathematical  model  is  reduced  to  a  computer  model  and 
computed  results  are  displayed.      The  design  of  a  pulse -position  modu- 
lator is  shown,    and  experimental  measurements  of  power  spectra  are 
compared  with  the  theoretically  determined  spectra.      The  theoretical 
determination  of  power  spectral  density  and  the  actual  measurement  of 
spectra  of  finite-time  samples  are  discussed. 

The  major  portion  of  the  research  reported  on  was  conducted  by 
the  writer  during  his  Industrial  Tour  at  Sylvania's  Electronic  Defense 
Laboratories  (EDL)  in  Mountain  View,    California. 
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1.  Introduction 

The  objective  of  this  report  is  to  relate  a  theoretical  approach 
for  the  determination  of  power  spectral  density,    to   show  computations  of 
power  spectral  densities  performed  for  a  particular  class  of  signals,    and 
to  compare  experimental  measurements  of  the  power  spectral  densities 
to  that  predicted  by  theory. 

The  power  spectral  density  of  a  pulse -position  modulated  (PPM) 
waveform  typically  consists  of  a  continuous  component  and  discrete  fre- 
quency components.      The  relative  amount  of  the  total  power  contained  in 
each  of  the  component  types  is  determined  by  the  modulation  process. 
The  distribution  of  power  in  the  spectrum  is  determined  by  the  modula- 
tion process,    the  pulse  shape,    and  the  average  pulse  repetition  rate. 

In  this  study,    a  mathematical  model  is  developed  which  determines 
the  power  spectral  density  of  PPM  waveforms.      The  mathematical  model 
is  then  reduced  to  a  digital  computer  program.      The  computer  program 
permits  the  computation  of  the  power  spectral  density  of  a  PPM  wave- 
form that  is  modulated  by  a  uniform  process,    a  Gaussian  process,   or  an 
exponential  process.      The  modulation  process  that  results  in  a  statistical 
correlation  between  pulses  within  the  pulse  train  with  respect  to  their 
displacements  (in  time)  is  developed  for  the  Gaussian  process  only.     Sta- 
tistical independence  between  pulses  is  assumed  for  the  other  processes 
considered. 

The  pulse  shape  is  assumed  to  be  rectangular  for  the  computer 
model.      The  pulse  width  and  the  average  pulse  repetition  rate  are  variable 


parameters  of  the  program.      The  pulse  amplitude  is  also  a  variable  param- 
eter,   thereby  allowing  computations  to  be  performed  for  combined  pulse- 
position  and  pulse -amplitude  modulated  waveforms.      By  setting  the  pulse 
position  modulation  equal  to  zero,    the  computation  of  the  power  spectral 
density  of  a  pulse  waveform  that  is  only  pulse -amplitude  modulated  may 
also  be  accomplished  with  the  computer  model. 

A  discussion  of  the  use  of  finite -time  interval  measurements, 
using  a  periodogram  technique  to  give  estimates  of  the  power  spectral 
density  of  a  process,    is  included  in  this  report.      The  periodogram  tech- 
nique is  used  to  determine  the  power  spectral  density  of  a  PPM  waveform. 
Since  the  characteristics  of  the  experimental  PPM  waveform  are  also 
used  as  an  input  to  the  computer  model,    a  comparison  of  the  theoretical 
and  experimental  results  is  included. 


2.  A  Discussion  of  Spectral  Analysis  (General) 

The  purpose  of  a  communication  system  is  to  convey  information 
in  some  form  from  one  place  to  another.      The  system  would  normally 
include  a  transmitter,    a  transmission  medium,    and  a  receiver  which 
recreates  a  reasonable  representation  of  the  original  message.      The 
message  creating  the  transmitted  signal  will  most  likely  be  one  of  a 
class  of  stochastic  phenomena;  that  is,   the  message  is  a  random  func- 
tion of  time.     Also,    the  noise  accompanying  the  signal  (introduced  by 
the  communication  system)  will  normally  be  a  stationary  Gaussian  ran- 
dom process.      The  analysis  of  such  signals  is  difficult  unless  one  re- 
sorts to  the  use  of  the  statistical  properties  of  the  signals.      The  power 
spectral  density  furnishes  us  such  an  analysis  technique,    as  it  is  rela- 
ted to  the  statistical  parameters  of  the  signal.     In  particular,    the  power 
spectral  density  is  related  to  the  covariance  (or  autocorrelation)  function. 

In  the  development  of  spectral  analysis  techniques,   one  normally 
considers  the  random  signal  to  be  a  stationary  process.     A  stationary 
process  is  one  in  which  the  statistical  parameters  of  the  signal  are  the 
same  everywhere,    independent  of  which  long  interval  over  which  they 
are  measured.     A  similar  concept  is  to  consider  a  process  as  an  en- 
semble of  signals  (an  ensemble  is  a  set  of  similar  things).     If  the  statis- 
tical averages  of  the  ensemble  measured  at  one  instant  of  time  are  the 
same  for  any  other  instant  of  time,   then  the  process  is  also  said  to  be 
stationary  in  time. 

If  a  single  process  of  the  ensemble  of  signals  has  statistics  over 


a  long  period  of  time  that  are  the  same  as  the  statistics  of  the  ensemble 
at  a  given  time,    then  the  process  is  said  to  be  "ergodic".     An  ergodic 
process  will  always  be  stationary,    but  a  stationary  process  does  not 
necessarily  imply  ergodicity  [6]  . 

In  this  paper  all  processes  will  be  assumed  stationary  unless 
otherwise  stated. 

To  arrive  at  the  statistical  parameter  used  to  obtain  the  power 
spectral  density,    consider  the  random  ergodic  process    y(t).      A  typical 
member  of  the  ensemble  is    y   (t).      The  ensemble  averaged  correlation 
function  (covariance  function  K(T))is  defined  as  the  second  joint  moment, 


or 


K(T)   =  EWt)}' 


where    E   denotes  the  expected  value  taken  over  the  ensemble.        The  auto 
correlation  function   R   (t)   of  a  typical  member  of  the  ensemble  is 


_n.    .         lim  1       P  n.   ,   n.     ,       .   , 

R   (T)   =  T-co   2T     J_T   y   "^   (t+T)dt 

In  general, 

K(t)    =   E{Rn(T)}  , 

where  the    E  j  R    (tH   here  is  the  autocorrelation  function  averaged  over  the 

The  statistical  definition  of  covariance  would  require  that  the 
square  of  the  mean   (iiy    )  be  subtracted  from  the  above  expression.      It  was 
decided,    however,    to  include  the  mean  in  this  investigation,    thus  permit- 
ting a  simplified  development  that  will  handle  signals  that  possess  a  non- 
zero mean. 


entire  (infinite)  ensemble.      This  average  should  normally  be  performed 
before  the  time  average  of  R(T). 

However,    for  the  ergodic  case 


Rn(T)    =   K(T)        (probability  1)  . 


The  mathematics  relating  the  power  spectral  density  to  the  autocorrela- 
tion (or  covariance)  function  were  first  developed  by  Wiener  and 
Khintchine.      The  Wiener-Khintchine  Theorem  states  that  the  power 
spectral  density   S(f)    is  the  Fourier  transform  of  the  covariance  func- 
tion,  or 

pOO 

S(f)    =    2  K(T)e"jUJ   dT,     f  £  o 

J-co 

Note  that  in  the  above  expression   S(f)    is  defined  for  positive  frequencies 

only,    hence  the  factor  of  two  in  the  equation.      This  is  obviously  the  more 

practical  case,    and  is  permissable  since   S(f)   is  an  even  function  (as  is 

R(t)).     One  can  obtain  the  autocorrelation  function  by  taking  the  inverse 

Fourier  transform  of  the  power  spectral  density  (again  assuming  ergo- 

dicity).      Thus 

i  r°° 

R(T)    =   -  S(f)cosuuTdf  (probability  1) 


o 


and 


p  uu 
S(f)    =    4  R(T)CosujTdT  (probability  1) 

Jo 

It  should  be  mentioned  that  while  a  signal  has  a  unique  power  spectral 


density,    the  converse  is  not  true.      All  phase  information  in  the  signal 
is  lost  in  S(f),    thus  the  original  signal  could  not  in  general  be  recov- 
ered from  its  power  spectral  density.      Furthermore,    a  particular   S(f) 
theoretically  represents  the  power  spectral  density  for  a  multitude  of 
signals. 


8 


3.  Finite -Time  Interval  Spectral  Analysis 

The  developments  of  the  previous  section  followed  the  typical 
theoretical  approach  to  signal  analysis  in  that  infinite  time  intervals 
were  assumed.     In  practical  analysis  situations,    such  intervals  are  of 
course  always  finite.      For  some  processes,    a  finite  investigation  in- 
terval can  be  considered  adequate  if  the  interval  is  large  with  respect 
to  the  fluctuation  time  of  the  process  under  investigation.     Since  finite 
time  measurements  of  spectra  of  PPM  waveforms  will  be  compared 
with  the  results  from  the  mathematical  model,    a  discussion  of  possible 
sources  of  error  due  to  analysis  of  the  finite-time  interval  observation 
of  the  signal  is  included. 

The  problem  of  finite  time  interval  analysis  of  a  stationary  ran- 
dom process  can  be  stated  quite  simply  as  follows.     A  sample  function 
of  the  random  process    y(t)    is  observed  during  an  interval   t   =    -T/2   to 
t   =    T/2.     Will  this  observation  contain  enough  of  the  characteristics  of 
the  process  to  give  a  reasonable  estimate  of  the  power  spectral  density 
of  the  process?    Note  that  the  random  process  under  consideration  is 
stationary.      The  stationarity  of  the  signal  becomes  more  meaningful 
(and  necessary)  when  considering  finite-time  analysis. 

To  develop  the  theory  for  finite-time  analysis,    again  consider 

the  ensemble    y(t),    except  now  let   y(t)    vanish  outside  the  interval    (-T/2. 

T/2).     An  amplitude  spectrum  for  a  member  of  the  ensemble  is 

poo  .  pT/2 

,rn...  n.   .        -juut  ,  n.   .    -juut  , 

Y    (f)       =  y   (t)    e   J     dt    =  y   (t)e   J     dt 

1         J-oo  X  J-T/2 


Now  the  finite -time  power  spectral  density  of  the  ensemble  is 

S(f)T    =   E-f!  Yn(f)TYn(f)T*  }    =    |r  E  {|Yn(f)T|2}  s 

where         denotes  the  complex  conjugate.      The  true  power  spectral  den- 
sity  S(f)    is  defined  in  terms  of  the  ensemble  average  over  an  infinite 
time  interval  (and  is  equivalent  to  the  definition  on  page  5).     Specifically: 

lim       „,„,  lim 


■w-£„««T«£„£*{|Y-(fLl2} 


1 
where  the  ensemble  average  is  carried  out  before  the  limit  is  taken. 

Again,    consider  the  sample  function  of  a  single  random  process. 
The  finite  time  transform  is  equivalent  to    Y    (f)       of  the  preceding  para- 
graph,   and  will  now  be  used  as  an  approximation  to  the  Fourier  trans- 
form of  a  single  stationary  process.      A  finite  power  spectral  density  of 
the  single  random  signal  may  now  be  defined. 

Sn(f,T)   =|   |Yn(f)T|2 

2,  3 
S(f,  T)    is  called  a  "periodogram",  and  has  received  considerable 

attention  due  to  the  fact  that  it  is  one  type  of  practical  measurement 

that  can  be  made  on  real  processes  (signals).     Observe  that  the  power 
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sity,  August  1962,   30-33 

3_ 

Grenander,    U.    and  Rosenblatt,    M.      Statistical  Analysis  of 

Stationary  Time  Series.     John  Wiley  &  Sons,    1957,    91-94 


spectral  density  for  the  time  limited  ensemble  is  the  expected  value  of 
the  periodogram  over  the  entire  ensemble,   or 


S(f)T    =   E{sn(fsT)} 


and 


«f>  -  £       «{s»«.T)} 


Thus,   one  measure  of  the  reliability  of  the  periodogram  is  the  fact  that 
as    T-»co,   the  mean  value  of  the  periodogram  does  equal  the  power 
spectral  density.     However,   as  shown  by  Abramson  and  Farison,    for 

some  very  common  random  processes,    the  variance  of  S(f,  T)  does  not 

2 
go  to  zero  as    T -»  oo.        The  fact  that  the  periodogram  (being  a  random 

variable)  has  the  desired  mean  value  but  may  vary  extensively  about  the 
mean  for  certain  processes  has  caused  some  skepticism  about  its  use 
as  a  spectral  density  estimator  [  5]  . 

It  has  been  shown  that  the  periodogram  is  not  a  consistent  esti- 
mate of  the  power  spectral  density  [  10]  .      However,    the  periodogram 
may  be  used  to  give  consistent  estimates  of  the  power  spectral  density 
at  specific  frequencies    f..     This  is  accomplished  by  the  use  of  weight- 
ing functions,   or  "spectral  windows"  about  the  frequency   f..     A  possible 
method  then  of  obtaining  a  consistent  estimator  of  the  power  spectral 
density  would  be  to  have  many  adjacent  windows,    and  at  the  end  of  the 
analysis  time    T   observe  the  response  of  all  windows  simultaneously. 
A  smoothing  between  the  point  estimates  observed  would  naturally  be 
desirable. 


Periodogram  calculators  fall  into  three  general  classes;  that  is, 
the  filtering,    computing  and  smoothing  may  be  performed  in  three  dis- 
tinctly different  ways.      These  are 

1.  Parallel  filter  banks, 

2.  Sweeping  analyzers, 

3.  The  coherent  memory  delay  integrator. 

The  parallel  filter  bank  technique  is,    as  the  name  implies,    a  bank 
of  narrow-band  filters,    each  with  a  common  input  (the  signal  to  be  ana- 
lyzed).     At  the  output  of  each  filter  is  an  energy  detector,    which  is 
periodically  sensed  by  a  sweep  circuit  that  is  connected  to  the  vertical 
deflection  plates  of  an  oscilloscope.      The  parallel  filter  bank  technique 
is  depicted  in  block  diagram  form  in  Figure  3-1. 


Input 


■> 


Spectral  window 
for   fj 


Spectral  window 


for   f 


N 


Energy 
Detector 


Energy 
Detector 


-> 


■> 


Sweeping 

Readout 

Device 


Output 


Figure  3-1.     Parallel  Filter  Bank' Periodogram  Calculator. 

The  number  of  filters  (N)  may  vary  with  the  frequency  range  (analyzer 
bandwidth)  and  the  resolution  required.      The  bandwidth  of  a  filter  is 
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approximately  the  inverse  of  the  rise  time.     Obviously,    the  bandwidth 
of  the  filters  determines  the  resolution  of  the  device.      The  "analysis 
time",    or  Period   T,   using  this  method  is  approximately  the  rise  time 
of  the  filter.      From  this,   one  may  observe  that  the  resolution  and  the 
analysis  time  are  essentially  the  inverse  of  each  other.      Consequently, 
if  one  designs  an  analyzer  for  fine  resolution,    then  the  signal  must  be 
observed  for  a  longer  period,    and  the  converse  is  also  true.      This  im- 
portant relationship  holds  for  periodogram  calculators  in  general. 

The  sweeping  analyzer  is  the  more  common  general  laboratory 
type  of  device.     With  this  device,   the  input  signal  is  mixed  with  a  sweep 
ing  local  oscillator,   the  output  of  which  is  filtered  through  a  single 
"spectral  window".     The  signal  is  then  detected  and  displayed  on  an  os- 
cilloscope.     A  block  diagram  for  the  sweeping  analyzer  is   shown  in 
Figure  3-2. 


Input 


Mixer 


Sweeping 
Local 
Oscillator 


Narrow 
Band 
Filter 


^      Detector 


■^    Output 


Figure- 3-2.     Sweeping  Analyzer  Periodogram  Calculator. 

A  disadvantage  of  this  type  of  analyzer  is  the  adverse  effect  of  the  tran- 
sient condition  at  the  input  to  the  filter  while  the  local  oscillator  is 
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sweeping.     If  the   sweep  rate  is  too  rapid,    the  filter  may  continue  to  ring 
due  to  this  transient  effect,    even  though  the  signal  is  no  longer  present 
at  the  input  to  the  filter,    and  the  resolution  is  no  longer  the  bandwidth 
of  the  filter. 

The  third  type  of  periodogram  calculator  to  be  discussed  is  the 
coherent  memory  delay  integrator.     This  device  is  somewhat  sensitive 
in  its  operation,    but  the  concept  is  unique  and  merits  consideration. 
The  block  diagram  of  the  device  is  shown  in  Figure  3-3. 


Input 


Readout  after 
N  loops 


Detector 


Delay 
T 


Unity  Gain  Storage  Loop 


Modulator 
(Phase 
Shifter) 


1 


Phase  Advance  Signal 
cos  [if  J  t 


-^»  Output 


Figure  3-3.     Coherent  Memory  Delay  Integrator 

The  local  oscillator  used  for  the  phase  advance  signal  must  be  precisely 
matched  to  the  delay  of  the  delay  line.  Furthermore,  the  loop  gain  must 
be  unity.     With  these  optimum  conditions  assumed  in  being,    the  response 

of  the  analyzer  to  a  single  frequency  is  a  main  lobes   with  side  lobe  pat~ 

2 

terns  very  much  like  a   (sin  x/x)      curve.      The  side  lobe  patterns  may 
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be  removed  in  some  instances  by  a  process  of  selective  "weighting"* 
in  which  the  input  signal  is  time  weighted  during  the  analysis  interval. 
As  an  example,    a  cosine  weighting  function  is  accomplished  by  ampli- 
tude modulating  the  input  signal. 

The  output  is  produced  by  a  linear  addition  of  a  large  number  of 
sinusoids  differing  from  each  other  in  frequency  by  small  increments. 
Phase  coincidence  of  all  signals  present  in  the  storage  loop  occurs  once 
(for  a  given  frequency)  during  the  observation  interval.      The  time  of 
this  phase  coincidence  is  linearly  related  to  the  input  frequency,   thus  a 
readout  of  the  response  after  the  proper  number  of  loop  recirculations 
yields  the  desired  periodogram.      The  one  operational  analyzer  of  this 
type  with  which  this  writer  is  familiar  is  known  under  the  trade  name 
"Simoramic". 
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4.  The  PPM  Spectral  Model 

Pulse -position  modulation  (PPM)  is  a  form  of  pulse  time  modu- 
lation in  which  the  modulating  process  causes  the  position  (in  time)  of 
the  pulse  to  vary  relative  to  its  unmodulated  time  of  occurrence.      The 
relation  of  the  pulse  to  the  modulation  process  is  normally  determined 
by  the  sampling  technique.      The  two  more  common  types  of  sampling 
are  "natural  sampling  and  "uniform"  sampling  [2]  .      For  natural   sam- 
pling,  the  occurrence  of  the  pulse  is  determined  by  the  instant  at  which 
the  leading  edge  of  a  sawtooth  sampling  wave  and  the  modulating  (mes- 
sage) voltage  are  equal. 

message  waveform  x([t) 

sawtooth  sampling  wave 
t 


resultant  PPM 
waveform   y£(t) 


Figure  4-1.     Generation  of  PPM  by  Natural  Sampling. 

When  uniform  sampling  is  used,    the  modulating  waveform  is 
sampled  at  equally  spaced  time  intervals.     Each  pulse  in  the  PPM  wave- 
form is  then  delayed  in  time  from  its  unmodulated  position  by  a  time    e 
which  is  proportional  to  the  sampled  voltage. 
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message  waveform  x{t) 
sampled  value  x^ 


resultant  PPM 
waveform   y1(y) 


Figure  4-2.      Generation  of  PPM  by  Uniform  Sampling. 

If  one  desires  to  relate  a  specific  modulating  process  to  the  re- 
sultant power  spectral  density  of  the  PPM  waveform,    it  is  necessary 
that  the  pulse  positioning  process    e   be  expressable  in  terms  of  the  mod- 
ulating process   x.      This  was  done  for  the  case  of  uniform  sampling; 
however,   the  attempt  to  do  this  for  the  natural  sampling  case  proved 
fruitless.      The  resulting  expression  for  the  natural  sampling  case  was 
in  the  form  of  an  unwieldy  series  of  expansion  which  could  not  be  reduced 
to  a  closed  form  nor  was  of  any  value  in  its  series  form.     Therefore, 
the  spectral  model  development  will  be  based  on  a  single  channel  PPM 
system  using  uniform  sampling. 

The  zero  order  hold  is  employed  by  the  uniform  sampling  tech- 
nique just  described.      The  model  can  be  adapted  for  use  with  other  sam- 
pling schemes.      For  a  Gaussian  process,   one  can  use  a  higher  order 
hold  in  the  model,      as  the  model  will  be  developed  to  handle  the  case     , 


An  example  would  be  a  technique  which  requires  that  the  value 
of   c,     be  linearly  related  to  the  average  of  the  previous   N    samples  of 
x(t). 
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for  which    e     ,e      are  correlated  Gaussian   bivariate  random  variables 
m      n 

for  a  finite   m-n.      The  correlation  for  this  hypothetical  case  can  be  in- 
troduced by  the  sampling  technique,    or  it  may  be  directly  related  to  the 
correlation  of  the  sampled  values  of  the  modulating  process,   or  both. 
Let  the  modulating  waveform  x(t)   vary  from  zero  to    A    volts. 
Define  the  maximum  allowable  deviation  of  the  pulse  from  its  unmodula- 
ted position  as    T    /L  seconds.      The  time  shift    e,     is  then 
o  k 

T  T 
o       .       .          __o 

Sk   "   LA  X(kTo)    "   LA  *k  ' 


The  pulse  time  shift    e      is  thus  linearly  related  to  the  sampled  amplitude 

of  the  message  waveform  x(kT    ).     The  variance  of   e,     is  also  linearly 

o  k 

related  to  the  variance  of  the  message  waveform  samples   x,  ,   or 

T  T  2 


\z  =  E{e~^e}2  =  E{Lix"d  ^}  >°: 


2        fTo\2      2 

\      -   Ila)   Qx     ■ 


The  standard  deviation  (or  RMS  deviation  from  the  mean)  of  the  pulse 

time  shift   e.     is 
k 

T 

a      =  -2,  a 
e        LA    x 


A  more  meaningful  expression  in  determining  what  the  actual 

RMS  deviation  is  relative  to  a  normalized  interval  will  be  defined  as    <J  . 

n 

a  a 

a     =  -*   =  -^ 
n        T  LA 

o 
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A  general  expression  for  the  power  spectral  density  of  a  PPM 
waveform  will  be  shown.  The  process  x{t)  producing  the  modulation 
will  be  considered  stationary.  The  development  of  the  equations  will 
follow  very  closely  that  of  Middleton  [8]  . 

A  pulse  function  U(t  »  nT     -  e)    shall  be  defined  as  having  any 
desired  shape  (with  maximum  value  of  one)  and  existing  only  for  a  dura- 
tion u  beginning  at  time    t  =  nT     +  e    .     Outside  of  this  specified  inter*- 

4 

val,   the  pulse  function  is  zero.     Define   y    "   as  the  peak  amplitude  of  the. 

i —  member  of  the  ensemble  of  PPM  waveforms  in  the  time  interval 

nT     +  e     <  t  <  nT     4 •  e     +  u.      A  member  of  the  truncated  ensemble 
on  on 

y    (t)  may  now  be  written: 

N 

*N  (t)    S     I    ynlu(t  "  nTo   "  0 

n=-N 

The  autocorrelation  function  of  the  PPM  waveform 


i.  .        lim  i.  . 

»  (t)   ■  N-co  %  (t) 


is    r\t). 


*V>  - 1£.  £  £'//nK>  ♦ T* 


N 
lim  1  r1  i      i 

N-»oo    (2N   +  1)T         L  ym  yn 

°m,n=-N 


I 


oo  .  , 

U(t  -  mT     -  e     X)U(t  +  t  -  nT     -  e     )dt 
o         m  on 

oo 


Now  let   X TT(f)   be  the  Fourier  transform  of  U(t). 
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pOO 

:u<f>  = 

-oo 


U(t)e"jUJtdt 


or 


riOO 

^-00 


t-a)e"jaJtdt    =   Xu(f)e"jaja 


Also 


oo 


J"* 


Xu(f)eJW"df 
•oo 


or 


oo 


Xu(f)ejU3(t"P)df 
-oo 


Letting   a  =  nT     +  £      -  T   and    P  =  mT     +  £        in  the  above  equations ,   we 
on  o  m 

have 


R(T)    = 


1 


lim 


N-oo  (2N  +  1)T 


S  .        .  poo 

n  =  -N  -OO 


o  m,n=-N 


riOO 


-OO 


Xu(f)ejUj(t"3)dfdt 


1 


lim 


N-oo  (2N  +  1)T 


N  .oo 

Y       y   y  x 

Li  n     m  J 

^-_M  "OO 


u 


(f)e 


-juup 


°  m,n=-N 


U(t  -  a)eJ     dtdf 


-oo 


L 


lim 
N-oo  [2N 


N 

I 

°  m,n=-N 


1 V  i      ifv    m    -J«* 

TTyT       L      WnJ     XT^(f)e 


'U 


Xu(f)*ejaJ0Cdf 
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R(t)  =— J        |x0(f)|   e 
o     -oo 


1 


litr. 


N-  oo     (2N  +   1) 


N 

I 

m,  n=-N 


.       .   Me 
yn  V  e 


-6     X)    juu(n-m)T  1 
n     m  o 

e  J 


df 


Now  is  it  necessary  to  determine  the  covariance  function  from  the  above. 
Recall  that  the  covariance  function  is  defined  as 


K(T)    =    e{rV)} 


Additionally,    recall  that  it  is  necessary  to  perform  the  ensemble  aver- 
age prior  to  taking  the  time  average  of  the  autocorrelation  function. 
Thus 


1 


oo 


|x|(f)|V*°T 

o      -co 


N 


1 


m,  n=-N 


Lists 


Nr+co    (2N  +   1) 


.    juu(e     -e       )   juu(n-m)T 
11  n       m 

y       y     e  e 

m     n 


°} 


df 


Now  observe  three  important  facts.      The  expected  value  of  a  finite  ser- 
ies it:  being  taken,    such  that    y  vanish  for       m    ,      n     >  N.     Since 

m,  n 

the  process  is   stationary,    the   statistical  relationship  of  any  two  random 

variables  displaced  the  same  distance   kT      in  time  is  the  same  for  all 

o 

such  pairs  of  random  variables.      Finally,    there  are  2N  +   1    such  like 
processes  in  the  double  series.      Define   k  =  m  -  n   in  the  series,    and  the 
covariance  function  may  be    x-cwritt  en 
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oo 


o      -oo 


lirh 
N-'oo 


N 

z 

k  =  -N 


E{yly2e 


juu(e    -Sj). 


juukT 


df 


t   -t  =kT 

2      1         o 


Now  that  the  mean  of  the  ensemble  has  been  taken,   one  may  apply  the 
limit,    and 


K(T)    = 


oo 


-oo 


oo 


T  lxu<f)l2    I 


k  =  -oo 


j.(c2-£l) 
Eiyly2e  J 


jtukT 


t    -t    =kT 
2      1         o 


e  df 


K(T)    is  now  in  the  form  of  an  inverse  Fourier  transform;  thus  the 
Wiener-Khintchine  Theorem  for  two  sided  density  functions  may  be  uti- 
lized at  this  point.      Let   S'(f)   be  the  two   sided  power  spectral  density. 


oo 
S'(f)    =   f     |Xu(f)|2    I      E{yiy. 


k=-oo 


juu^-ejh  juukTQ 


kT 


-00  <  f  <  CO 

Later,    the  power  spectral  density  shall  be  defined  for  positive  fre- 
quencies only  (one  sided),    and  since    S'(f)    is  an  even  function,    the  one 

sided  power  spectral  density  is    S(f)    =    [  2S'(f )]  u    (f),    where    u    (f)    is 

o  s 

the  unit  step  function.      The  unit  step  function  is  defined  as 

0  ,     f  <  0 

u    (f)    =    1/2  ,     f  =  0 
s 

1  ,     f  >  0 
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Define    S     ' (f )   as  the  two  sided  spectral  density  for  the  typical  sampling 


pulse,   or 


lxu(f)| 

V(f)    =    T '      "°°   ^f  ^°°  * 


S     '(f)    is  a  continuous  even  function  of  frequency.      The  one  sided  spec' 

tral  density  for  the  typical  sampling  pulse  would  then  be 

2|X    (f)|2 
Su(f)    =   2Su'(f)   =  - ,     f  *  0 


=    0 


,     f  <  0  . 


Thus,    the  power  spectral  density  may  be  rewritten 

juukT 


CO 


s«(f)  =  su«(f)   2,  E|yiy2e  }     e        ,-oo^f^oo. 

k=-oo  kT 

o 


The  spectral  model  to  be  used  assumes  that  the    y      are  indepen- 
dent from  sample  to   sample,    and  independent  from  the    £   process. 
Thus,    the  equation  that  is  used  for  the  generalized  spectral  model  is 


S'(f)    =    Su'(f) 


2       _2 

y    +  y 


CO 

I 

k  =  -co 


E    e 


juu(e    -gj) 


{.  ■  2  '  } 


juukT 


kT    =t    -t. 
o      2      1 


,       -OO    ^  f  ^  co 


Observe  that  the  model  can  also  be  used  for  pulse -amplitude  modulation 
(PAM)  in  which  the    y      represent  the  voltage  level  of  the    k —    sample  of 

.K. 

the  message  and  the  pulse  is  subjected  to  a  random  jitter  (in  time), 
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possibly  due  to  equipment  temperature  fluctuations  or  varying  transmis 
sion  path  length.      The   y      would  of  necessity  be  independent  random 
variables  for  the  generalized  spectral  model  just  developed.      The     £    , 
however,    as  will  be  shown  in  the  further  development  of  this  PPM 
model,    can  either  be  independent,   or  correlated  random  variables  with 
a  Gaussian  distribution. 
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5. 


A  Discussion  of  the  Statistics  in  the  Spectral  Model 


As   seen  in  the  development  of  the  spectral  model,    certain  sta- 
tistics of  the  two  processes    e  and  y,    together  with  the  pulse  shape,    in 
essence  determine  the  power  spectral  density.     Since  independence  was 
assumed,    the  statistics  of  the  pulse  amplitude    y      are  quite  simply  de- 
termined,    as  all  that  is  involved  is  the  mean  of  the  square  of  the  pro- 
cess (y    )  and  the  square  of  the  mean  (y     ).      The  statistics  of  the    £    pro- 
cess   are  not  as  simple,    especially  since  independence  of    £      will  not  be 
assumed  for  all  cases.      The  statistics  of  the    e    process  are  directly  re- 
lated to  the  message  process  x(t).     It  has  already  been  shown  that  the 
variance  and  the  standard  deviation  of  the  two  processes  are  linearly 
related.      Furthermore,   the  probability  density  functions  (pdf)  of  the  two 
processes  are  of  the   same  form  if  one  considers  continuous  distribu- 
tions only  (i.e.,    if  x    is  Gaussian,    e    is  also  Gaussian).     Specificallyj 
if  p       (x)    is  the    pdf  for  x,    then 


ekU)  = 


o  pXl  w  J 


k  v    e 


To  show  that  this  is  true,    recall  that    E  =  (0    /a    )x  .     Now 

£      x 


P[e     <  e]     =   P 
k 


a   £ 
x 


k       a 


But 


a  £ 

k        a 


a  e/a 

x  £ 

-co 


p      (x)dx 
Xk 


Since  the   pdf  is  the  derivative  of  the  distribution  function,    then 
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/  1      d 

k 


«a  e/o 

J  Px    (x)dX 

-co  k 


from  which  one  can  obtain  the  desired  result 


a  /a  e 

P.   (£)   =  7T  P, 


e.  a     x  x.  \    a 

k  e        k  \     e, 

To  avoid  carrying  unnecessary  terms  in  equations  already  complex  with 
terminology,    a  time  axis  shift  of  the  PPM  waveform  will  be  assumed 
such  that  the  mean  of  the  process    e    is  zero.      That  this  simplification 
may  be  performed  with  no  effect  on  the  resultant  power  spectra  is  in- 
tuitively obvious,    but  nevertheless  has  been  previously  verified  by  the 
writer  in  each  case  considered. 

To  evaluate  the   statistics  of  the  process    £    that  appear  in  the 
equation  for  the  power  spectral  density,    one  must  determine  the  joint 
characteristic  function   G(juu.,    juu    )  of  a  bivariate  process.      The  joint 
characteristic  function  is  defined  as 

G(jwr   ju>2)    =   E|e  |    -   E|e  j    . 

This  joint  characteristic  function  will  be  evaluated  for  the  Gaussian 
bivariate  case.      For  the  remainder  of  the  cases  included  in  the  compu- 
ter model,    statistical  independence  from  interval  to  interval  of  the    e 
process  was  assumed.      For  the  independent  case,    the  problem  reduces 
to  the  determination  of  the  square  of  the  absolute  value  of  the  charac- 
teristic function  of   e. 
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E^e 


juu(e2-e1); 


E    e 


juue  \  I  2 


G  (jut)) 

Cm 


The  equation  for  the  power  spectral  density  now  reduces  to 


S'(f)    =SuI(f) 


2  +  -2  \r 
y     +  y      |G 


co 


»l2  I 


jwkT, 


k  =  -co 

k^O 


CO    ^   f    ^  CO 


But 


CO  CO  / 

-co  o      -co      >  o 


Recall  that   STT(f)    =    2S'(f)   for   f  ^  0   and  S(f)    =    2S'(f)u    (f).      Then 
U  s 


Su(f)<y2  + y2  |G£((ju)|Z 


CO 


1  + 


r    I 

o     k  =  0 


6    f  - 


u   (f) 

s 


for  the  case  in  which    £     and  y,     are  independent  from  sampling  interval 
to  sampling  interval.      As  mentioned  previously,    only  positive  frequen- 
cies are  considered. 
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5A.  The  Message  as  a  Gaussian  Process 

The  purpose  of  this  subsection  (and  the  following  subsections)  is 

to  fit  a  specific  process  into  the  generalized  spectral  model.     The  pro- 
- 
cess considered  may  either  be  the  message  process    x   or  the    e   process, 

as  the  similarity  of  the  two  processes  has  been  previously  shown.      For 
the  Gaussian  case,    the  development  will  allow  for  a  statistical  correla- 
tion to  exist  between  I   adjacent  samples  of  the  process.     More  speci- 
fically,   this  development  will  allow  for  the  spectral  model  to  handle 
Gaussian  processes  for  which  the  autocorrelation  function   R(T)   exists 

for   x(t)   up  to     I T I    =   IT    ,    and   R(t)    =    0   for    I T I  >  IT   .      This  is  practi- 

o  o 

cal  because  a  purely  random  process  yields  an  autocorrelation  function 
that  approaches  zero  as    T   increases  [1], 
The  expected  value 
jMe2-€j) 


-{•        2     l  } 

kT 


o 


is  the  joint  characteristic  function  of   e    ,  e^    for  a  given  kT      -  t      -  t., 

The  joint  characteristic  function    G(ju).,juu    )  as  previously  defined  fits 

the  case  now  under  development  when   uu     =  -<■*}_.     It  was  decided  to 

carry  out  the  general  development  of   G(juu    ,  juj    ),    substituting  into  the 

final  result  ou  =  uu     =  -*uu   .      Thus, 

{  ju>(e    -e   )  j(di  e  +w  e   ) 

E{e         2  }    =   E{e       l    '      2    Z}t    GO^.jO,) 


The  expected  value  is  determined  as 
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oo     poo    (jijUjej+juu    e2) 
G(juu1,juu2)  J        e  p(eis  e2)de1de2 


-oo       -oo 


The  bivariate  Gaussian  pdf  is 

h  [-1/20   (1-p  )]  (gj    -2pe1e2+e22: 

^•i-^  = :  2  a    2  e 

2ttQ    V1   "  P 
where 


u,      =    0    and  a        =   a        =    o 
e  £1  £2 


For  simplicity,    substitute 


£1    =X  £2    =    * 


2a2   =   a  Vl   ~  P2   =  b 


Then 


2         2 

iCo    juu  jMi|t    /ab    ) 


1J~2'         -rrab 


-oo 


e-(l/ab   )(X    -Z^H    1    dxdt 


-OO 


But 


-(l/ab2)(\2-2p\i(r)   jV 


oo 

oo 

X 


=    eP 


2t2/ab2p°°     -(l/ab2)RX-p*)2]     jmr 

e  e  qa. 


-oo 
which  is  in  the  form  of  a  constant  multiplying  the  characteristic  function 

2  2 

of  a  single  Gaussian  process    \,   with   mean   |X    =  pt  and  cr        =  ab    /2. 

A.  A. 

The  characteristic  function  of  a  single  Gaussian  process  may  easily  be 
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determined  (with  the  aid  of  tables  of  integrals),    and  the  last  integral  is 


simply 


2    2         2 
i —       p   |i    /ab 
bva-rr    e 


2      2 
j pij/cu ,  -(ab   fju  1    /4) 


The  joint  characteristic  function  then  becomes 


rl.    .    ,       i      -(abV2/^    r°° 

-co 


(p2-l)/(abV   j{'VP'Ul^ 


di|j  I 


Define 


OU 


3    =    m2  +   ptDj 


1  /       ab 


2  V       2 

P      "  ] 


and  this  integral  also  takes  the  form  of  the  single  characteristic  function 
of  the  process    i|j ,    with  the  mean   u     =  0 

1 


/     u2       %. 

it-       -(ab  (ju1  ;/4 
G(J(Wl'JaJ2)  =  QWl    6 


rj      -(ab   (JUj    /4) 


V2~V 


2     2 

oo     ~(l/2a      )i|r      jua^ 

e  e  diji 


-oo 


•«^V/2)' 


By  re  substituting  the  values  for    0\  ,    a,    b,    and  uu-  one  obtains 


G(j(JU1<iuu2)    =    e 


2  2  2 

-(a      /2)(uu,    +2p(w1oj_+fju_    ) 


Recall  that  for  this  development,   uu   =   U)      =    -UK.      Thus 


G(j(ju1juj    )    =    G(fju)    =    e 


2    2/i       x 
-a     uu   (1-p) 

e 


Applying  the  above  to  the  spectral  model  for  a  given  value  of    kT    ,    we 


28 


have 


E     e 


juu(e    -Sj) 


-a  2/[l-p(kT   )] 
e  o 


kT    =t    -t. 
o     2     1 


That  the  expression  is  a  function  of  the  autocorrelation  of  the  process  e 
(or  x)  may  be  determined  by  observing  that  p  is  the  correlation  coeffi- 
cient for  a  given  kT      in  the  development.      Considering  all   kT    ,    p(kT    ) 

o  o  o 

is  the  normalized  autocorrelation  function. 


p(T)    = 


R(T) 

R(0) 


Now  let   R(t)  =  0   for   k  >  I.      The  power  spectral  density  for  this 


case  is  then 


f 


S-(f)    -    Su'(f)/y2  +  y"2 


~*        -a£2oj2[l-p(k)]     jU)kTQ 

L  e 

k=-I 

k/0 

V/"?1)j*To  V  jajkTo 

+    e  )        e  +       )        e 


k=(I+l) 


k  =  -co 


But    p(kT    )  =  p(-kT    ),    as    R(t)  is  an  even  function,    and 
o  o 


S'(f)    =   Sy'U)  l  y2  +  y2 


•    I    /    -acV[l-P(k)]     -c/\ 

Z    )    I  e  -  e  I    cos  uukT 

Lk=P  ' 


oo 


jiukT 


2    2, 
-<J      fJU    / 

+    e  f   i   1  +      ;        e 

>  k="-oo  ^ 


and  the  spectral  model  for  the  Gaussian  case  (for    f  ^  0)  becomes 


Su(£!     y2  +  y2 


I    /   -a£2ju2[l-P(k)]         -o*^ 

)     (  e  -  e  J  co  s  'jokT 


2  L  l 
.  k=i 


+    e 


2    2 

■a     a> 
e 


°  k=o  ° 


u   (f) 
s 
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5B.  The  Message  as  a  Uniform  Process 

The   specific  spectral  model  will  be  developed  for  the  case  in 
which  the  message  determining  the  modulation  of  the  PPM  waveform  is 
a  uniform  process.      Statistical  independence  of  both  the    e    and  the    y 
processes  will  be  assumed.      For  this  assumption,    the  spectral  density- 
was   shown  to  be 


S(f)    =    STT(f){y2  +   y2   |G    (UU)|2 


CO 


1  +  f  I6  £-f )  k(£> 


o    k  =  0 


The  problem  is  therefore   reduced  to  that  of  determining  the  character- 
istic function   G   (uu)   of  a  uniform  process. 


G  (uu)   = 
e 


r>  to     . 

eJ      p      (e)de 
-co  k 


As   shown  previously,    if  x    is  uniformly  distributed,    then    e    is  also  uni- 
formly distributed.      For  the  case  of  uniform  distribution  of  x,    the  total 

possible  deviation    e(max)  allowed  is    T    /L.    seconds  (or  ±(T    /2L)  from 

o  o 

the  zero  mean  assumed).      Therefore, 


k  o 


2L  2L 


Thus, 


or 


G   (uu) 
e ' 


T    /2L 

o  L      juue  , 

~r  e       de 

~[V  /2L)         o 
o 


sin 


G    (uu) 
e 


uuT 
o 

2L 


ijuT 
o 

2L 
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The  power  spectral  density  for  the  uniform  case  then  becomes 


uuT 
.  sin  — 
S(f)    =   S    (f){/+yZ' 


uuT 
o 

~2l7 


1  + 


00       / 

n     1 r>        \ 


O    k  =  0      N  O- 


U     (f) 

s 
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5C.  The  Message  as  an  Exponential  Process 

The   specific  spectral  model  will  be  developed  for  the  case  in 
which  the  message  determining  the  modulation  of  the  PPM  waveform  is 
an  exponential  process.     Of  interest  for  this  case  is  that  the  speech 
process  is  approximately  exponentially  distributed,    especially  in  the 
region  of  medium  to  large  amplitudes  [1]  . 

For  the  exponential  case,    statistical  independence  will  again  be 
assumed  for  both  the    £  and  y   processes.      As  was  the  case  then  in  sec- 
tion 5B,    the  problem  is  reduced  to  determining  the  characteristic 
function  of  the  process    e. 

The  exponential  probability  density  function  is 

j  -(V2/a   )|e| 

p      (e)    =— —   e 

from  which  one  can  determine 


\1      =    E(e)    =    0 
e 


Var(e)    =   E(e2)  -  |i  2    =   E(e2)    =   a   2 

c  e 


Define 

V2 


a  = 
a 
e 


Then 


,   .       a     -a  e 
p     (e)  =  -  e 

£k  2 
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The  characteristic  function  is 


CO 

G   (uu)  ejUJ£p      (e)de 

-co  k 


a 

2 


e(j(ju+a) 
e  de  + 


-co 


CO 


e(juu-a) 
e  de 


or 

ufl    +  a 

The  power  spectral  density  of  the  PPM  waveform  for  the  exponential 
case  is  then 


(uu    +  a  ) 


-  1  + 


fH-f 


u    (f) 
s 
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6.  The  Computer  Model 

The   PPM  spectral  model  was  developed  in  section  four,    and  the 
model  was  adapted  to  specific  modulation  processes  in  section  five. 
The  resultant  equations  are  quite  formidable  if  manual  calculation  is  re- 
quired to  determine  the  power  spectral  density  for  various  cases.      For- 
tunately,  the  tedium  of  manual  calculations  of  such  proportions  is  not 
necessary  nowadays,    as  the  digital  computer  handles  such  tasks  with 
great  facility.      A  computer  model  will  now  be  developed  that  will  incor- 
po  rate  all  cases  previously  discussed,    plus  the  unmodulated  case. 

The  computer  model  developed  relates  the  pulse  positioning  pro- 
cess   e    in  terms  of  the  original*modulating  process.      This  allows  the 
user  to  have  a  more  direct  insight  as  to  the  effects  of  the  modulating 
process  on  the  power  spectral  density  of  the  PPM  waveform.     Most  of 
the  parameters  of  the  model  are  variable  at  the  discretion  of  the  user. 
An  exception  is  the  number  of  points  of  continuous  spectra  calculated 
between  each  possible  discrete  frequency  spectrum.      Since  five  points 
in  each  such  interval  are  adequate  for  a  graphic  plot,    this  quantity  has 
been  fixed  at  five  to  avoid  possible  confusion  by  future  users  of  the 
model. 

It  should  be  noted  here  that  the  following  development  of  the  com- 
puter model  may  not  at  times  appear  to  be  the  simplest  approach.      The 
approach  taken  generally  yields  simplification  of  computer  Fortran  oper- 
ations to  be  performed,    such  as  eliminating  the  calling  of  a  library 
function  for  performing  an  exponentation  when  a  simple  multiplying 

34 


operation  will   suffice.     It  is  admitted,    however,    that  a  few  of  the  "sim- 
plifying" techniques  employed  did  not  necessarily  make  the  program 
more  efficient.      The  development  will  include  these  also,    as  they  have 
become  a  part  of  the  final  computer  model. 

While  no  extreme  efforts  were  applied  toward  making  the  com- 
puter program  efficient,    this  aspect  was  certainly  kept  in  mind  while 
the  program  was  being  developed. 

The  program  was  written  in  Fortran  IV  language  for  the  IBM 
7  0  40  at  Sylvania's  Electronic  Defense  Laboratories  (EDL).     With  some 
obvious  card  changes  (for  each  library  function  used),    the  program  can 
be  adapted  to  the  Fortran  language  of  the  CDC   1604  at  the  Naval  Post- 
graduate School. 

The  total  run  time  (including  card  reading  and  program  printout) 
for  computing  four  separate  cases  in  one  program  run  from  a  binary 
deck  was  timed  at  two  minutes  on  the   7040. 

From  the  previous  developments,    it  was  decided  to  base  the 
computer  model  on  the  mathematical  model  expressed  in  the  following 
form  for  frequencies  greater  than  zero.      The  case  for    d-c  (zero  fre- 
quency) will  be  developed  separately. 

S(f)    =    Su(f){y2+  y2[Sum(f)   -    |G£(uu)|2]} 

CO 

+  su(f,{?2|GeW|2f  Hf"f)   } 

°k=l  ° 

The  first  half  of  the  expression  is  the  continuous   spectrum  and  the 
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second  half  is  the  spectrum  concentrated  at  discrete  frequencies.      The 
term   Suml'i:    is  defined  as 


I     .    -(a£w)2[l~p(k)]  ~(o^)Z 

2  )    f  e  -  e  I  cos  kuuT 

k=l 


I 

a(f)  2  £   ( 

1,-1  V 


The  term   Sum(f)    is  included  in  the  model   so  that  the  Gaussian  correla- 
ted process  may  be  handled.      Note  that  when   Sum(f)    equals  zero,    the 
model  reduces  to  that  derived  for  the  case  in  which  all  processes  were 
assumed  independent.      The  term    G    (uu)    is  simply  the  characteristic 
function  of  the  independent  random  process    e,    and  fits  as  shown  into 
the  chosen  representation  even  when  the  Gaussian  correlated  case  is 
considered.      That  this  is  true  may  be  determined  from  the  fact  that  the 
characteristic  function  for  the  independent  Gaussian  case  (mean  zero) 
is  the  square   root  of  that  determined  for  the  Gaussian  bivariate  case 
when    rho  (p)    is  zero,    or 

~[(aeuo)2/2] 

G   (uu)    =   e 

e 

(for  Gaus  sian  process  with    p    =    0). 

All  terms  in  the  model  have  now  been  adequately  defined  except 
for  the  power  spectral  density  of  the  "typical  sampling  pulse".      Recall 
in  a  previous  development  that  this  pulse  could  be  of  any  shape  with  a 
maximum  amplitude  of  one.      Also  recall  that  this  amplitude  of  one  does 
not  restrict  the  amplitude  of  the  pulses  in  the  PPM  waveform,    as  the 
PPM  pulse  amplitudes  are  determined  by  the    y   process.      Thus,    the 
"typical  sampling  pulse"  is  simply  a  description  of  the  shape  of  the 
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pulses  of  the   PPM  pulse  train.      For  this  computer  model,    a  rectangular 
pulse  will  be  assumed,    the  pulse  being  of  duration   u   and  of  unit  ampli- 
tude. 


-u-J 


Figure  6-1.      The  Typical  Sampling  Pulse  for  the  Computer  Model. 


As  previously  shown, 


Su(f)  -  ~    |xu(f)|2. 

o 


Here 


oo 
Xu(f)  U(t)e"JUJtdt    = 

-co 


from  which  one  may  determine 


urn 

2  /  Sln  ~r 

ST,<f>  V 

U  uuu 

2 


u/2 

■(u/2) 


d)e-jUJtdt. 


To  simplify  the  computer  program  (so  that  some  of  the  time  consuming 
operations  are  minimized),    the  above  equation  may  be  rewritten  as 


c     .  .         4u        (1   °  cos  uuu) 
Su(f)    =   —  2        " 

o  (uuu) 


The  discrete  frequencies  at  which   S(f)    is  to  be  computed  are 


n 


n         5T 


or   oo 


2Trn 
n         5T 


0  ^  n  ^  200 
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To  calculate   S(0),   one  must  apply  the  value  of  the  unit  step 
function  (one -half)  at  zero  frequency  (the    d-c    term),     Substituting  uu  -   0 
into  the  equations  developed  thus  fars   one  obtains  for  all  cases 
STT(0) 


where 


STT(0)  2 

U  u 

2  T 

o 

The  sequence  of  computation  of  the  power  spectral  density  may 
be  determined  from  the  flow  chart  (Figure  6-4). 

The  computer  model  will  determine  two  hundred  values  of  S(f) 
up  to  and  including    f  =  40/T   .      This  includes  all  40  possible  values  of 
S(f)  concentrated  at  discrete  frequencies  (delta  functions)  within  the  fre- 
quency range  considered. 

A  list  of  the  significant  symbols  used  in  the  program  and  their 
relationship  to  the  terms  in  the  mathematical  model  is  shown  in  Table 
6-  1.  The  computer  program  itself  with  one  page  of  output  is  included 
as  Figures  6-2  and  6~3. 
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Table  6  =  1.      Computer  Program  Symbol  Identification. 


Symbol 


Equivalent  Symbol  in  Mathematical  Development  or 
Symbol  Definition 


A 

ALFA2 

CFSQ(N) 

CONT(N) 

DELFN 

FL 

F(N) 

I 

J 

NR 


a     (used  lor  exponential  case) 

|G   (uu)|2 
e 

Value  of  continuous  spectrum  at    f  =  N/5T 

o 

Value  of  discrete  spectrum  component 
L 

Frequency;  F{N)  =  N/5'T      cps 

I    (upper  summation  index  for  Gaussian  case) 

None;  this  index  specifies  modulation  process.     See 
comments  in  Figure  6-2. 

Number  of  separate  cases  to  be  computed  this  program 
run 


RHO{K) 

P(k) 

SIGMA 

a 

X 

S(N) 

S(f); 

SU(N) 

S   (f) 

u 

SUM(N) 

Sumi 

T 

T 
o 

TEMP(K) 

Part 

U 

U 

YBSQ 

2 

y 

YSQB 

y 

S(f);  Power  spectral  density 


Summation  for  Gaussian  correlated  case 


Partial  value  of  summation  of  SUM(N) 
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C  PROGRAM  PPM 

C  THIS  PROGRAM  CALCULATES  THE  POWER  SPECTRAL  DENSITY  OF  A  PPM 

C  WAVEFORM  THAT  HAS  AS  ITS  MODULATING  SOURCE  ONE  OF  THE  FOLLOWING 

C  EITHER  AN  INDEPENDENT  GAUSSIAN  PROCESS  OR  A  GAUSSIAN  PROCESS 

C  WITH  UP  TO  TEN  CORRELATION  VALUEStFOR  TAU»T,2T , . .., 1QT) ,  AN 

C  INDEPENDENT  UNIFORM  PROCESS.  OR  AN  INDEPENDENT  EXPONENTIAL  PROCESS 

C  J=l  FOR  GAUSSIAN, *2  FOR  UNIFORM, =3(UNM0DULATED) ,=4  FOR  EXPONINTIAL 

c    Format  8  statement  allows  for  reading- record  in "following  order 

C      U,T,YBSQ,YSQB,SIGMA,A,FL, J, I, AND  10  VALUES  OF  RHO 

D I MENSI ON  RHO (TO) , SUM  ( 2  00 ) , C  FS  Q ( 200 ) , SU ( 200  J , ST200 S ,  T EMP ( 10 )   , 

1F(200),    C0NT(200)  

6   FORMAT    <4H   F=   F10.3,    6H   S(F)=   E12.4,    9H  CONTIF I =~ET2  .  47 

7  FORMAT    (/5X,1HU,IIX,IHT,9X,4HYBSQ.8X,4HYSQB,6X,5HSIGMA, 6X, 1HA, 10X, 
12 HF L , A X ,2 9H    2ND    ROW    IS    J , I ,       RHO ( K= 1 , 10 ) / 1 

8  FORMAT    (7F11. 6/(213, 10F6. 31) 

9  FORMA T    ( 4H F»    F TO.  3 I  a r  6H \   S  ( F  J F*    E 12 .4 ) 

PI=    3.14159265 

NR    =    4 

300  PRINT    7 
REAO"8ru;frYB"SQ"tYSQB,SIG"MAVAfFL,J,r«CRHOIKr,"K»lVl01 
PRINT    8,    U,T,    YBSQ,YSQB,SIGMA,A,FL,J,I,(RH0(K),K=1,10) 
DC I    30irN=r,200 

301  SUM(N)    =    0.0 

GO    TO    (1,2, 3,4), J 

1    C=    (0.4«PI«SIGMA/(FL»A))««2 

DO    11    N =1 ,200 

FL02    =    N 

C  C : =C •  FL  02  •  F  L 02 

IF    (CC.LT.23.0)    GO    TO    12 

CFSQ(N)    =    0,0 

GO    TO    11 
12   CFSQCnT    «    1.0/ (EXP  (CO) 
11    CONTINUE 

IF    ( f  .EQ.    <5)    GO    TO  TOO 

DO    18    N«l,200 
FLUJ   »  U 

TEM    »    0.0 


DO  15   K»T,I 

FL04    =    K 

CU   »   C*FL03«FLO3«(1.0-RH0(K)l 

IF    (CD. LT. 23.0)    GO    TO    14 

TERPTkl  «  0TT5 

GO   TO    15 

14  TEMPTKJ    *    (T.O/(EXP(CD))-CFSQ(N)  J  •(C0S(0.4»PI  »FL03«FL04  )  ) 

15  TEM    »    TEM*TEMP(K) 
1 8    SUM (NT  »    2 70 •TEM 

GO    TO    100  

~TTU  ^~0  7Z»  PT7FT 
DO   25    N=l,200 
FL06»N 
PUN=PU«FL06 
25    CFSQ(N)«((SIN(PUN))/PUN)«»    2 

GO    TO    100      

3    DO    35    N=l,200 
35    CFSQ(N)    «    1.0 


Figure  6-2.  Source  Program  for  PROGRAM    PPM,    Page  1 
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GO  TO  100 
A  ALFA2  =  2.0*(FL«A/(T*SIGMA) )»»2 

PE  =  (0.4«PI/T)*»  2 

DO  45  N=l,200 

FL05  =  N 
45  CFSQ(N)=(ALFA2/(PE*FL05»FL05  +ALFA2))«»  2 
100  AA  =  4.0*U*U/T 

B=  0.4»PI»U/T 

DF=  0.2/T 

DO  110  N=i,200 

FLOl  =N 

F(N)=FLOl  ♦  OF 

G  =  B»FL01 
110  SU(N)  =  AA*(1.0-C0S(G) )/(G»G> 

SUZERO  =  AA/4.0 

SZERO  =  SUZEROMYSQB  -YBSQ*( 1.0-1 .0/T  )  ) 

FZERO=  0.0 

DO  120  N=l,200 
120  SIN)  =  SU(N)«(YSQB*YBSQ*(SUM(N)-CFSQ{N) ) ) 

DO  125  N=5,200,5 

CONT(N)  =  SIN) 

DELFN  =  SU(N)«YBSQ«CFSQ(N)»1.0/T 
125  S(N)  =  CONT(N)*DELFN 

PRINT  9,  FZERO,  SZERO 


DO  210  K=l,196,5 
K3=K+3 

200 
210 

DO  200  N=K,K3 

PRINT  9,  F(N),S(N) 

N=K+4 

PRINT  6,F(N),S(N),C0NT(N) 

NR  =  NR--1 

IF  (NR.GT.O)  GO  TO  300 

C 

THIS  RUN  FOR  10  KC  SAMPLING 
STOP 

END 

U  =  0.1T 


Figure  6-2,   Page  2 
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START 


NR  =  #  of  vA'aveforms  considered  in  this  run] 


(Read  waveform  "\ 
characteristics  y 


waveform  \ 
.c teristics  J 


rint  wavefo  rm  N 

charai 


Set   SUM(N')  -   0 
N  -  1,  200 


J  -  4 
(exponential) 


J  =  1 
(Gaussian) 


J  =  2 
(uniform) 


Compute 
CFSQ(N) 
N  -  1,200 


CFSQ(N)  =  1.  0 
N  =■!,  200 


Compute 
CFSQ(N) 

N '  =  1,  200 


Compute    SU(N) 
N  =  1, 200 


Compute 
CFSQ(N) 

N  =  1.  200 


No 


Yes 


Compute 
SUM(N) 
N  =  1,  200 


Compute  SU(O) 
Compute  S(O) 


Compute  continuous  spectrum,    CONT(N 
N  -  1,200 


i 


Compute  discrete  spectrum,  DEL.FN 

for    F  ■=  N/5T 


Compute  total  power  spectral  density 
S(N)  =  CONT(N)  plus  DELFN 


Print    F(N),   S(N):  also  print  CONT(N) 
when   S(N)  contains  discrete  spectrum 


No 


?NR  =  NR  -   1 


Figure  6-4.      Flow  Chart  for  PROGRAM  PPM. 
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7.  Computed  Power  Spectral  Densities,    Results 

This  section  contains  the  displays  of  the  power  spectral  densities 
computed  for  certain  combinations  of  pulse  widths,    modulation  processes, 
and  varying  degrees  of  pulse  deviation  (or  RMS  deviation  from  the  mean 
position).      The  legend  accompanying  each  graph  indicates  the  pertinent 
data  regarding  the  process  for  which  that  power  spectral  density  was 
computed. 

The  first  graph  (Figure  7-4)  is  the  power  spectrum  of  an  unmo- 
dulated pulse  train.  Observe  that  with  no  modulation  of  the  pulse  train, 
there  is  no  continuous  component  in  the  power  spectrum. 

Three  different  normalized  autocorrelation  functions  were  used 
for  the  Gaussian  case.      For  the  process  which  has  very  little  correlation 
between  samples,    the  approximation  used  was  the  linear  autocorrelation 
function  indicated  by  Figure  7-1. 


->T 


Figure  7-1.     Limited  Correlation  Between  Samples, 


An  approximation  to  the  autocorrelation  function  of  a  process  with  a  high 
degree  of  correlation  existing  for  several  adjacent  samples  is  shown  in 

Figure   7-2. 

14 


1  IT 


>    T 


Figure, 7-2.     Samples  Highly  Correlated. 

An  approximation  to  the  autocorrelation  existing  for  band-limited  Gaussian 
noise  (0~W  cps)  is   shown  in  Figure  7-3.     In  Figure  7-3,    the  Gaussian 
noise  is  assumed  to  have  been  sampled  at  ten  times  the  highest  frequency 
present;    or  at  10W  samples  per  second. 


>    T 


Figure  7-3.     Approximate  Autocorrelation  for Oversampled 
Band-limited  Gaussian  Noise. 

Two  examples  are  shown  for  which  the  pulse  amplitudes  are  not 
held  constant  (Rayleigh  distributed),    in  addition  to  being  pulse  position 
modulated.     It  may  be  observed  that  in  these  cases  the  continuous   spec- 
trum begins  with  an  appreciable  value  at  zero  frequency,    which  is  a 
characteristic  of  pulse  amplitude  modulated  waveforms. 
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To  allow  for  an  adequate  display  of  the  continuous   spectrum,    a 
scale  factor   was  chosen  that  does  not  permit  the  showing  of  the  discrete 
power  spectrum  to  scale.      The  values  of  the  discrete  power  spectral 
lines  are  numerically  indicated  for  each  line  that  exceeds  the   scale  of 
the  graph. 
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8„  Spectral  Analysis  of  Actual  PPM  Waveforms 

A  comparison  of  the  results  from  the  spectral  model  was  made 
with  the  results  determined  experimentally.      The  objective  of  this  aspect 
of  the  investigation  was  to  determine  if  short  time  interval  analysis 
would  yield  a  reasonable  estimate  of  the  power  spectral  density  of  PPM 
waveforms. 

In  order  to  generate  PPM  waveforms  comparable  to  those  asso- 
ciated with  the  spectral  model,    a  pulse  position  modulator  was  designed 
and  built.      The  circuit  and  performance  characteristics  of  the  modula- 
tor are  included  in  Appendix  A.      The  unmodulated  pulse  train  generated 
has  a  pulse  repetition  rate  of  one  kilocycle,    and  a  pulse  width  of  195 
microseconds. 

The  spectrum  analyzer  (periodogram  calculator)  used  for  the 
experimental  measurements  is  known  by  the  trade  name  "Spectran". 
The  Spectran  is  one  of  the  parallel  filter  bank  types  of  analyzers  des- 
cribed in  section  3.     It  contains  480  parallel  filters,    and  has  a  resolution 
of  25  cycles  per  second  (cps)0      The  analysis  time  is   18  milliseconds. 
The  first  two  frequency  bands  of  the  Spectran  were  used  for  the  PPM 
analysis      covering  20  to  3852  and  3852  to  7684  cps,    respectively.      The 
use  of  two  bands  accounts  for  what  appears  to  be  a  slight  discontinuity  in 
the  pictures  to   be  presented  later  in  this  section. 

The  Spectran  was  calibrated  by  inserting  a  sine  wave  input  of  an 
arbitrary  frequency  (two  kilocycles)  and  recording  the  analyzer  response 
for  measured  power  inputs.      The  input  power  was  considered  the  square 
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of  the  RMS  voltage  input,    which  is  equivalent  to  the  spectral  model  defi- 
nition 01  powtr  into  a  one  ohm  resistor.      The   resulting  calibration  curve 
is  included  as  Figure  8-1.     Unfortunately,    the  calibration  curve  is  only- 
accurate  for  spectra  at  discrete  frequencies  (delta  functions)  as  a  conse- 
quence of  having  been  based  upon  data  taken  for  discrete  frequency  in- 
puts.     The  calibration  curve  will  not  yield  an  accurate  value  of  the  con- 
tinuous spectrum  at  a  specific  frequency   f  .      Each  of  the  480  "spectral 
windows"  of  the  Spectran  has  a  bandwidth  of  about  10  cps.      Therefore, 
when  a  continuous   spectrum  is  measured,    the  measurement  obtained  at 
a  frequency   f      is  actually  the  proportionate  power  of  all  frequencies  that 
are  passed  through  the   spectral  window  that  passes  the  frequency    f   . 
Further,    the  noise  level  of  the  measuring  equipment  prohibited  accurate 
measurements  of  the  continuous  portion  of  the   spectrum.      Consequently, 
for  the   experimental  results,    a  determination  of  the  exact  values  of  the 
continuous  spectrum  must  be  attempted  by  a  different  procedure.      The 
general   shape  of  the   continuous  spectrum  was  determined  by  the  Spectran 
analysis,    however     and  is  displayed  in  Figure   8-3.      A  discussion  of  the 
background  related  to  the  photographs  in  Figure  8-3  is  considered  neces- 
sary in  order  that  they  may  be  properly  interpreted.      As  predicted  by 
the  computer  model,    the  value  of  the  power  spectral  density  at  any  fre- 
quency of  the  continuous  spectrum  is  several  orders  of  magnitude  less 
than  that  of  the  larger  discrete  components.      In  addition,    the  predicted 
average  levels  of  the  continuous  spectrum  is  about  the  same  as  the  noise 
level  output  of  the  measuring  equipment  (the  Spectran).      Thus,    if  one 
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visually  subtracts  the  noise  level  and  the  effects  of  the  discrete  compo- 
nent as  determined  from  the  top  picture  in  Figure  8-3,    then  a  fair  esti- 
mate of  the  shape  of  the  continuous   spectrum  may  be  gleaned  from  the 
bottom  picture. 
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Figure  8-1.      Spectran  Calibration  Curve 

2 
A  plot  of   (V_wo)      at  a  discrete  frequency  vs  Spectran  out' 
RMS 

put  response.     (Spectran  input  attenuator  set  at  28  db.  ) 
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In  order  to  measure  the  value  of  the  continuous  power  spectrum 
at  discrete  frequencies,    a  wave  analyzer  with  a  tunable  6  cycle  bandwidth 
was  employed.        The  PPM  waveform,    modulated  by  a  random  noise 
(assumed  Gaussian),    was  used  as  the  input  to  the  wave  analyzer.      True 

RMS  voltage  measurements  of  the  analyzer  response  were  made  at  100 

2 
cycle  increments.         The   block  diagram  of  the  experimental  setup  is 

shown  in  Figure  8-4. 


1  kc 

Pulse  Train 


G.R.  Type  1390A 

Noise  Gen. 

2  0  kc  Bandwidth 


PPM 

Modulator 


Narrow  Band 
Tunable  Filter 


True  RMS 
Voltmeter 


HP-302A 


Fluke  9 10- A 


Figure  8-4.      Technique  for  Experimental  Measurement  of 
Continuous  Spectrum. 

Assuming  an  ideal  rectangular  response  of  unit  amplitude  for  the 
filter,    the  total  power  to  the  RMS  voltmeter  is 

f   4-3 


Power 


Or 


S(f)df 


T    -3 
o 


where    f      is  the  center  frequency  of  the  filter.     It  shall  be  assumed  that 

The  analyzer  is  the  HP-302A,    which  has  an  equivalent  noise 
bandwidth  of  6  cycles  with  the  center  frequency  adjustable  throughout  the 
audio  frequency  range. 

2 

True  RMS  voltage  measurements  were  made  with  a  "Hi-damped" 

Fluke  Model  910-A  voltmeter,    which  is  a  power  sensitive  instrument 
calibrated  in  volts  RMS. 
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there  are  no  extreme  changes  in  the  power  spectral  density  in  any  6  cycle 
interval     r  ru  j  ■    S(f)  may  be  approximated  by  a  constant   S(f   ).      The  meas- 
ured power  is  then  approximately 


Measured  Power  =  6S(f    )  . 

o 


S(f    )  may  now  be  determined  from  the  square  of  the  true  RMS  voltage 

measurement,    which  is  the  measured  power  (to  a  one  ohm  resistor), 

2 

V 

s(£  )  =  -5MS_ 

o '  6 

A  plot  of  the  power  spectral  density  measured  in  the  manner  just 
described  is  shown  in  Figure  8=5.      The  few  inconsistent  measurements 

recorded  are  attributed  to  this  writer's  inability  to  "time  average"  a 

3 

slightly  fluctuating  voltmeter  needle  by  visual  observations.        For  com- 
parison purposes,    the  theoretical  plot  of  power  spectral  density  is 
included  as  the  dashed  line  in  Figure  8-5. 

The  measurement  of  the  discrete  spectrum  was  a  simpler  matter, 
as  output  levels  well  above  the  equipment  noise  levels  were  being  meas- 
ured.     As  previously  mentioned,    the  calibration  chart  applies  to  the 
measured  results  of  the  discrete  spectra,    from  which  may  be  determined 
a  reasonably  accurate  measured  power  level.      Unfortunately,    the   smooth- 
ing of  the  Spectran  "readout"  is  not  perfect,    resulting  in  the  occurrence 


3 

For  a  discussion  of  the  variance  of  the  output  of  a  square  law 

device  (the  same  type  as  the  measurement  fluctuations  encountered),    see 
page  2  56,    Davenport,    W.    B.      and  W.    L.    Root,    An  Introduction  to  the 
Theory  of  Random  Signals  and  Noise,    McGraw  Hill,    19  58. 
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of  peaks  and  valleys  when  one  sweeps  across  the  spectrum  with  a  con- 
stant levt  .     .  pu1    signal.      Thus,    even  the  discrete  frequency  power  meas- 
urem<  nts  w  11   have  some  equipment  measuring  error  introduced. 

The  discrete  spectrum  displays  are  shown  in  Figure  8-2.      From 
F  gure   3-2  and  theoretical  computations,    an  experimental  versus  theore- 
tical    comparison  was  made  and  is  included  in  Table  8-1.      The  PPM  de- 
vice is  that  described  in  Appendix  A,,    modulated  with  a  wide  band  Gaussian 
noise   source. 


Discrete  Power  Spectral  Density 


Dis crete 
Frequency  j 


Unmodulated  Pulse  Train 


Modulated  PPM 

2       ~~2~ 
a     =  0.  112;  y     =  y     =  2.  72 

n 


Experimental        _,  . 

,,-,  ^.  „  J  Theoretical 

(From  Fig.  8  =  2) 


Experimental 
(From  Fig.  8-2) 


Theoretical 


Table   8-1.     Experimental  and  Theoretical  Values  of 
Discrete  Power  Spectrum. 


74 


9  5 

*     ° 


h 

e 

W  ' 

0  >s 

o  d 

«  u 

6  J 


09 


V 

a 

•0 

ID 

K. 

ti 

ti 

U 

3 

rt 

n) 

* 

* 

c 

vO 

-* 

CO 

J, 

CO 

00 

io 

rj 

•— » 

o 

V 
A 
H 


■ 

JJ 

tt     a 

flj 

k        m 

* 

*    * 

CO 

«M           sO 

«r> 

^        oa 

•                  • 

o 

o        o 

2 

0 

-t-> 

* 

u 

«) 

X 

o 

o 

a,     U 

* 

T* 

3  ■* 

0 

Q 

^     O   ifl 

•a 

V     =) 

-a 

4-> 

Pow 

ntin 
to  5. 

o 

o 

o 

H 

H 

U 

wor  -* 


9.  Discussion  of  the  Results  and  Conclusions 

A  reasonable  question  to  be  answered  now  is  "What  specific  sta- 
tistics of  the  modulation  process  contribute  to  the  shape  of  the  power 
spectral  density?"     This  is  one  question  that  is  answered  in  this   section. 
Other  items  discussed  in  this  section  are: 

1„       The  "general    shape"  of  the  power  spectral  density  of  PPM 
waveforms. 

2.  The  power  spectral  density  as  a  function  of  the  modulating 
process. 

3.  The  parameters  that  determine  the  relative  amount  of  the 
total  power  that  is  contained  in  the  discrete  and  continuous 
spectra. 

4„       The  effect  of  correlation  between  pulses  in  the  pulse  train 
with  respect  to  their  relative  positions  in  time. 

5.       Some  general  comments  on  the  experimental  determination 
of  power  spectral  density. 

The  shape  of  the  power  spectral  densities  of  PPM  waveforms 

r      .2 

appears  to  be  generally  of  the  form  (sin  x/xj     if  one  considers  only  the 
frequencies  beyond  the  first  lobe.         The  portion  beyond  the  first  lobe  has 
been  observed  to   be  mostly  continuous  spectrum  for  the  modulation  pro- 
cesses considered.      The  reason  for  this  can  be  determined  from  the  equa- 
tion upon  which  the  computer  model  is  based.      The  characteristic  function 

1 

will  be  assumed  that  all  pulse  waveforms  discussed  in  this 

section  consist  of  rectangular  pulses  except  when  otherwise  stated. 
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for  all  modulation  processes  considered  "damps  out"  toward  zero  value 
as  the  frequ*     >  y     s   increased,,    resulting  In  a  power  spectral  density  at 
higher  frequencies  that  is  essentially  the   same  as  the  spectrum   S    (f)  of 
the  typical   sampling  pulse. 

The  frequency  at  which  the  first  lobe  is  terminated  (where  the 
envelope  of  power  spectral  density  first  becomes   zero)  is  the  reciprocal 
of  the  pulse  width.      The   shape  of  the  first  lobe  appears  to  be  quite  depen- 
dent upon  the  modulating  process.      Figures  7  =  5,    7-8,    and  7-14  are 
representative  of  the  spectrum  of  PPM  waveforms  that  differ  in  the  modu- 
lating process;  I.  e.  s    a  uniform  process,    an  exponential  process,    and  a 
Gaussian  process,,    respectively.      The  biggest  difference  in  these  power 
spectral   densities  is  In  the  distribution  of  power  in  the  discrete  frequency 
components.        A  close  look  at  the  Figures  7-5,    7-8,    and  7-14  addition- 
ally reveals  a  distinct  difference  in  the  continuous  spectrum  within  the 
first  lobe.      As  was  the  case  with  the  discrete   spectrum,    this  difference 
may  be  attributed  to  the  type  of  statistical  process  used  to  modulate  the 
PPM  waveform. 

The  question  asked  at  the  beginning  of  this  section  will  now  be 
answered.      For  the  cases  in  which  statistical  independence  between 
samples  is  assumed,    the  square  of  the  absolute  value  of  the  character- 
istic function  is  the  only  statistic  of  the  modulating  process  that  influ- 
ences the  power  spectral  density  of  the  PPM  waveform.      For  all   cases 

2 
The  frequencies  at  wh'ch  the  discrete  components  exist  are 

multiples  of  the  reciprocal  of  the   sampling  rate,    or  the  average  pulse 

repetition  rate. 
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considered  in  this  investigation,    the  absolute  value  of  the  characteristic 
function  is   smply  the  Fourier  transform  of  the  probability  density 
function  of  a  process  with  zero  mean,      The  process  with  a  non-zero 
mean  yields  a  characteristic  function  that  is  equal  in  magnitude  to  the 
characteristic  function  of  the  process  with  a  zero  mean,    but  with  a  phase 
difference.      Taking  the  absolute  value  removes  this  phase  difference, 
thus  yielding  the  same  results  as  for  the  process  with  zero  mean. 

For  the  case  in  which  statistical  independence  between  samples 
is  not  assumed,    the  power  spectral  density  is  influenced  by  the  joint 
characteristic  function  of  the  modulating  process.      The  effects  of  corre- 
lated samples  will  be  discussed  in  greater  detail  in  the  latter  part  of  this 
section. 

The   relative  amount  of  the  total  power  in  the  continuous  spectrum 
of  a  PPM  waveform  with  pulses  of  fixed  duration  and  amplitude  appears 
to  be  determined  by  the  degree  of  fluctuation  of  the  pulse  shift.      This  ob- 
servation may  be  confirmed  by  comparing  figures  based  upon  like  modu- 
lating processes  but  with  different  standard  deviations  (or  different  total 
pulse  shift  permitted  in  the  uniformly  distributed  case).      Comparing 

Figures   7-11  and  7-14,    for  which   a      is  0.  05  and  0.  1,    respectively,    one 

n 

may  observe  that  an  additional  21.  6  percent  of  the  total  power  is  included 
in  the  continuous  spectrum  for  the  greater  fluctuating  pulse   shift.      For 
the  uniform  case,    the  PPM  fluctuations  may  be  increased  to  the  maximum 
level  by  permitting  the  maximum  pulse  excursion  (delay)  to  be  equal  to 
the  full  sampling  interval    T      less  the  pulse  width   u.      This  corresponds 
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to    L  -  T    /{T      -  u)    in  the  PPM  models  developed.     When  this  is  done, 
o         o 

the  discrete  spectrum  becomes  quite  small,    and  eventually  disappears 
for  very  narrow  pulse  widths.      The  latter  observation  may  be  confirmed 
by  direct  calculations  with  the  mathematical  model;  the  former  observa- 
tion may  be  confirmed  from  Figure   7-5.     In  Figure  7-5,    the  maximum 
pulse  excursion  within  a  sampling  interval  is  permitted.      The  pulse  width 
is  0„  IT    ,    and  the  aggregate  power  in  the  discrete  components  is  reduced 
to  less  than  one  percent  of  the  total  power  (excluding  the    d-c    component). 
Another  factor  which  influences  the  relative  amount  of  total  power 
contained  in  the  continuous  spectrum  is  the  pulse  width.      For  a  given 
modulation  process,    as  the  pulse  width  is  decreased,    the  percentage  of 
the  total  power  contained  in  the  continuous  spectrum  increases.     In  con- 
trast to  the  considerations  of  the  preceding  paragraph,    in  which  the  total 
power  of  the  PPM  waveform  is  constant  (independent  of  the  fluctuations), 
if  the  pulse  width  is  decreased  the  total  power  decreases  unless  a  corre- 
sponding increase  in  pulse  amplitude  is  made.      A  plot  of  the  percent  of 
total  power  in  each  component  as  a  function  of  the  pulse  width  is  shown 
in  Figure  9-1  for  a  single  modulation  process  (the  independent  Gaussian 

case  with   CT      -  0.  0  5).      Data  for  this  plot  was  obtained  from  the  computer 

n 

model. 

The  effect  of  statistical  correlation  between  sample  values  of  the 
modulation  process  tends  to  center  the  continuous  spectrum  about  the 
discrete  components.     It  is  interesting  to  note  that  the  relative  amount  of 
the  total  power  in  the  discrete  components   remains  essentially  the  same 
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for  the  correlated  and  the  uncorrelated  cases  considered,    yet  the   shape 
of  the  continuous  spectrum  is  radically  changed  (see  Figures  7-11,    12, 
and  13s    for  example).      However,    in  the  limit  as  the  correlation  function 
of  the  modulation  process  approaches  unity,    the  PPM  waveform  is  essen- 
tially modulated  by  a    d-c    level,    and  the  power  spectral  density  is  all  at 
discrete  frequencies.      In  fact,    the   spectrum  becomes  that  of  an  unmodu- 
lated pulse  train  (Figure   7-4). 

The  correlated  case  is  a  practical  consideration  for  typical  PPM 
pulse  trains  encountered.      For  most  PPM  systems,    there  is  time  divi- 
sion multiplexing  of  several  channels  of  data  pulses,    which  are  further 
interleaved  with  synchronizing  (sync)  pulses.      The  sync  pulses  are  trans- 
mitted at  a  periodic  rate.      The  exact  periodicity  of  the   sync  pulses  re- 
sults in  a  contribution  to  the  discrete  components  of  the  power  spectral 
density.      Further,    the  data  pulses  are  normally  statistically  related  to 
the  sync  pulses.      The  mean  position  of  the  data  pulse  is  related  to  the 
sync  pulse  by  some  fixed  interval  for  each  process.      The  degree  of  cor- 
relation between  data  and  sync  pulses  is  set  by  the  data  pulse  variance 
about  its  mean  displacement.      Finally,    adjacent  data  pulses  may  be 
correlateds    as  was  the  case  for  the  Gaussian  correlated  process  of  the 
simpler  single  channel  PPM  model  used  for  this  investigation. 

The  verification  of  the  theoretical   results  by  taking  experimental 
measurements  was  discussed  in  section    8.      The  measurement  of  the 
continuous  spectrum  of  the  power  spectral  density  can  be  quite  arduous. 
The  difficulty  is  partly  due  to  the  fluctuations  of  the  modulating  process 
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(or  the  degree  of  modulation)  varying  in  a  probabilistic  manner  during 
the  measuring  intervale      This  effect,    of  course,    can  be  minimized  by 
averaging  the  measurements  over  a  longer  measurement  interval.      Fur- 
thermore,   the  noise  introduced  by  complex  measuring  equipment  can 
obscure  the  low-level   signals  encountered  in  the  continuous   spectrum. 

The   slight  discrepancy  between  the  measured  and  the  calculated 
power  spectral  densities  displayed  in  Figure  8-5  may  be  largely  attributed 
to  the  non-linearity  of  the  pulse-position  modulator  for  the  larger  modu- 
lating voltages.      Consequently,,    the  standard  deviation  of  the  pulse  shift 
e    is  not  as  large  as  that  determined  from  the  measured  standard  devia- 
tion of  the  modulating  process.      This  results  in  a  smaller  continuous 
component  and  corresponding  increases  in  the  discrete  components  as 
compared  to  the  calculated  power  spectral  density,    which  was  based  upon 
a  linear  modulation  scheme,, 
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APPENDIX  A 
The  Pulse -Position  Modulator 

A  pulse-position  modulator  was  designed  and  constructed  so  that 
experimental  measurements  could  be  made  on  the  type  of  PPM  waveform 
used  for  the  mathematical  model.      The  design  specifications  were  to 
build  a  device  that  sampled  a  modulating  process  at  a  periodic  rate  and 
generated  a  pulse  train,    the  pulses  of  which  being  delayed  in  time  by  an 
amount  linearly  related  to  the  sampled  value.      To  accomplish  this,    it 
was  decided  to  build  a  sampler  which  had  as  its  "hold"  circuit  an  RC  net- 
work with  a  time  constant  that  resulted  in  an  approximately  linear  decay 
during  the  interval  under  consideration.      This  linearly  decaying  voltage 
was  then  compared  to  a  "zero  modulation"  reference.     When  the  two 
voltages  were  equal,    a  pulse  was  transmitted  to  a  monostable  multivi- 
brator.     The  monostable  multivibrator  was  used  to  give  a  pulse  output  of 
the  desired  shape,    and  to  allow  easy  variation  of  pulse  width  and  pulse 
amplitude.      See  Figure  A-2  for  pictures  taken  of  the  unmodulated  and  the 
modulated  pulse  trains  at  the  multivibrator  output. 

The  circuit  diagram  is  shown  In  Figure  A-l„      The  variable  parairi' 
eters  of  the  circuit  are  the  pulse  width  (change    Rj,    the  pulse  height 
(change    V      ),    and  the  sampling  rate  (change    R,    and  pulse  rate  of  the 
pulse  generator).      The  sampling  rate  is,    of  course,    the  average  pulse 
repetition  rate  of  the  PPM  waveform. 
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Horiz:     0.  2  msec/cm 
Vert:       1  V/cm 


Unmodulated  Pulse  Train 


Horiz:     0.  2  msec /cm 
Vert:       1  V/cm 


CT 


0.  08 


Pulse  train  modulated  with  wideband  Gaussian  noise. 


Figure  A-2.     Performance  of  the  Pulse -Position  Modulator. 
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